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Abstract 

We consider the following partial integro-differential equation (Allen-Cahn equation 
with memory) 

rt 

'0 

where e is a small parameter, h is a constant, f(4>) is minus the derivative of a double 
well potential and the kernel a is a piecewise continuous, differentiable at the origin, 



e 2 4> t = l\{t-t') [e 2 A«^ + /(0) +eh](t')dt'. 
Jo 
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scalar- valued function on (0,oo). The prototype kernels are exponentially decreasing 
functions of time and they reduce the integrodifferential equation to a hyperbolic one, 
the damped Klein-Gordon equation. By means of a formal asymptotic analysis we show 
that to the leading order and under suitable assumptions on the kernels, the integro- 
differential equation behave like a hyperbolic partial differential equation obtained by 
considering prototype kernels: the evolution of fronts is governed by the extended, 
damped Born-Infeld equation. 

We also apply our method to a system of partial integro-differential equations which 
generalize the classical phase field equations with a non-conserved order parameter and 
describe the process of phase transitions where memory effects are present, 

u t + e 2 <p t = Jo oi(t - t') Au(t')dt' 

< 

j e 2 <k = So a2(t - f) [ e 2 A<p + f(<j>) + e u }(t')dt', 
where e is a small parameter. In this case the functions u and 4> represent the temper- 
ature field and order parameter respectively. The kernels a\ and 02 are assumed to be 
similar to a. For the phase field equations with memory we obtain the same result as 
for the generalized Klein-Gordon equation or Allen-Calm equation with memory. 



1 Introduction 

In this manuscript we consider the following partial integro-differential equation 

e 2 t = f ait - t') [e 2 A0 + fU) + e h](t') dt', (1) 
Jo 

where <f>(x,t) is a field, called an order parameter, defined in Q x [0,T], where Q C R 2 and 
T > 0, with Dirichlet and Neumann boundary conditions. The parameter e is assumed to 
be small, e 1, f(<f>) is a real odd function with a positive maximum equal to 0*, a negative 
minimum equal to —<f)* and precisely three zeros in the closed interval [—b, b] located at 



and ±6, where 6 is a positive constant. For simplicity we will consider 6 = 1. The operator A 
is the Laplace operator. The kernel a is assumed to be a piecewise continuous, different iable 
at the origin and scalar-valued functions on (0, oo) satisfying additional conditions to be 
described later. 

Specific cases of equation (|J) are a(t) = S(t), the Allen-Cahn equation: 

e 2 i = e 2 A0 + /(0) +eh, (2) 

which is the simplest model of phase transition with a non-conserved order parameter, and 
the damped Klein-Gordon equation (used, e.g., in the theory of long Josephson junctions): 

e 2 <p tt + 7 e 2 (fH = e 2 A0 + /(</>) + e h, (3) 
which is obtained for an exponentially decreasing kernel of type [?] 

a(t) = e" 7 *. (4) 

Equation (^|) can be thought of as a phenomenological equation based on energetic penaliza- 
tion driving the evolution of the system toward equilibrium states. More specifically, let us 
call F e = Jr2[^€ 2 (V4>) 2 — V((p)]dx, with /(</>) + e h = dV/d<p the free energy of the system. 
The functional derivative of the free energy, 5F e (<p)/5(j) is considered as a generalized force 
indicative of the tendency of the free energy to decay towards a minimum. The equation 
is obtained by assuming that the response of <fi to the tendency of the free energy to decay 
towards a minimum is given by [?, ?], 

/■* SF 
e% = ~ a(t- t') °-± (f) dt'. 
Jo ocp 

A system of integro-differential equations closely related to (|1|), the phase field equations 
with memory, is [?, ?] 
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u t + e 2 4> t = Jo ai(t - £') Au(f)df 

(5) 

^ e 2 4 = / * a 2 (t - t'[ e 2 A0 + + e it ](*')<**' 

in f2x [0, t] where Q G R 2 and T > with Dirichlet and Neumann boundary conditions for the 
temperature field, u(x,t) and the order parameter, <j)(x,t), respectively and for u vanishing 
initially. The parameter e and /(</>) are as for ([I]) and the kernels a\ and a 2 are similar to a; 
i.e., piecewise continuous, differentiable at the origin, and scalar-valued functions on (0, oo). 

For ai(t) = a 2 (t) = 5(t), system (^) gives rise to the classical phase field equations with a 
non-conserved order parameter 

u t + e 2 (f)t = Aw, 

(6) 

e 2 t = e 2 A0 + /(0) + e M . 

System (|5|) also generalizes the hyperbolic phase field equations with a non-conserved order 
parameter: 

u tt + e 2 4> u + 7i ut + e 2 7i 4>t = a Am, 

(7) 

^ e 2 <p tt + e 2 7 2 t = e 2 A0 + /(</>) + e u, 
reducing to them when the kernels a\ and a 2 are exponentially decreasing functions; i.e., 
cii(t) = ai e~ 74 *, % = 1, 2, with and ji are non- negative constants. System (|?p is obtained 
by differentiating both equations in (|5|), rearranging terms and rescaling by means of the 
transformation t — > a 2 2 t, 72 — > a 2 72, 71 — > a 2 71 and a := ai/a 2 . 
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The phase field equations with memory (|5|) describe the process of phase transitions where 
memory effects are present. For a readable description of the classical phase field equations 
see [?] . The first equation in (|5|) is based on the balance of heat equation for a non- Fourier 
process in which the expression for the heat flux is given by a convolution in time between 
the temperature gradient and the kernel a x [?, ?]. The second equation is obtained in a 
similar way as equation ([[]) [?]. 

In what follows we assume that there exists a solution <f)(x,t) of ([]]) or {u(x, t), <p(x, t)} of 
(H) defined for all small e, every x G Vt and every t G [0,T] which contains an internal layer. 
We also assume, for such solutions, that for all small e > and all t G [0,T], the domain Q 
can be divided into two open regions f2+(i;e) and fi_(£, e) with a curve T(t;e), separating 
between them. This interface defined by 

r(t;e) := {x G Q : (f>(x,t;e) = 0} , (8) 

is assumed to be smooth, which implies that its curvature and its velocity are bounded 
independently of e. The function cf) is assumed to vary continuously across the interface, far 
from the interface tending to 1 when x G fi+(£; e), -1 when x G Q-(t, e), with rapid spatial 
variation close to the interface. The problem is to derive a closed equation for the evolution 
of the interface asymptotically valid as e — > 0. 

The latter problem has been studied formerly by means of formal asymptotic analysis for 
PDEs (|) and (|) [?, ?, ?], as well as for systems © and (0) [?, ?]. For @ and ©, fronts 
evolve according to the mean curvature flow equation [?, ?, ?], 

V = K, (9) 

where v is the normal velocity of the interface and k its curvature. If y = S(x, t) is the 
Cartesian description of the interface, equation (|9|) reads 



5 



For the undamped Klein-Gordon equation, (|3|) with 7 = h — 0, Neu [?] proved that the 
evolution of the interface is governed by the Born-Infeld equation 

(1 - S 2 t ) S xx + 2S x S t S xt - (1 + Si) S tt = 0. (10) 

For (|H) with 7 > and for (0) Rotstein et. al. [?, ?] showed that fronts move according to 
an extended version of the Born-Infeld equation given by 

(1 - a S 2 ) S xx + 2 a S x S t S xt - a (1 + Si) S tt - 



_ 7 St (1 + S 2 X - a S 2 t ) -v(l + S 2 x -a S 2 )* = 0, (11) 
where the parameter v (proportional to h in (|3|)) is defined later. 

In local (geometric) coordinates equation flTT|) reads 



- - 9 +-f v = k + v (l-av 2 )^. (12) 
1 — a V s 

where v, Vt and k are the normal velocity, normal acceleration and curvature of the interface 
respectively. Equation ([□]) or (|12|) have been studied in [?], and in [?] for 7 = v = . 

The case of integrodifferential equations ([]]) and @ is still less investigated. For (j^) and for 
kernels which are Laplace transforms of suitable functions satisfying 

roc roo roc 

I a(t) dt < 00, / a(s) ds < 00 and / a(s) s ds < 00 

Jo Jo Jo 

where a(s) is the inverse Laplace transform of a(t), Rotstein et al. [?, ?] showed that the 
equation governing the motion of the interface is also (|TT| ) where in this case 



4(0) 



7 /r,\2' 



In the present manuscript we show that in the asymptotic limit studied and for a certain class 
of kernels, which will be made clear below, equation (|l|) is equivalent to (^) from the point 
of view of fronts motion; i.e., the evolution of interfaces is described by the same equation, 
(0) with 

a = (^J a(r) dr^j , 
where «(r) is the inverse Laplace transform of a(t), and 

7 = a 2 ( [ a(r) r dr) . 



We also find the equivalence between (|5|) and (|7]) [?] for a more general class of kernels, 
described below. In Section 2 we describe the class of kernels considered in this work, give 
some examples and explain our strategy in dealing with them. In Section 3 we derive the 
equation of motion for (|1|). We present the derivation treating the equation in Cartesian 
coordinates and representing the kernel as a Laplace tranform of a suitable function. For 
kernels represented as a Fourier transform of suitable functions the derivation is similar. In 
our derivation we assume that the interface has no oscillations in the sense that there are no 
points on the interface for which the velocity vanishes; i.e., the front is an advancing front. 
This assumption is not necessary for the case of exponentially decreasing kernels; i.e., for (|3|). 
In Section 4 we derive the equation of motion for (|5|) for kernels which are inverse Fourier 
transforms of suitable functions. Here we treat the problem in local coordinates. 



2 Description of the class of kernels considered 

In order to deal with the memory integral in deriving the equations of front motion for either 
([[]) or (H) we represent the kernels a(t) and cii(t), i = 1,2, as Laplace transforms of suitable 
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functions: 



POO 

a(t) = / air) e~ T 1 dr, (13) 
Jo 

or 

poo 

aAt) = / aAr) e~ T 1 dr. (14) 
J o 

for i = 1,2. Substituting (0) into (0), rearranging terms and defining 

X(t;r)= /V'^ [e 2 A0 + /(0) + e (15) 

where it is understood that x is also a function of the space variable x and it is defined for 
all x G Q, we obtain 

e 2 </>t = Io° «( r ) X(^ r ) 

(16) 

Xt (t; r) + r X (t] r) = e 2 A0 + /(0) + e /i, 
for all r G [0, oo). Note that the second equation in fllBD satisfies x(0, r) = for all x E Q 
and for all r G [0, oo). We can see that, in spite of the fact that there is an integral involved 
in the first equation, that integral sums over the parameter r not involving the time variable 
t as in ([!]). For the phase field equations with memory, substituting ( |H| ) into ([]), rearranging 
terms and calling 

X(t;r)= fe^{t-t>) [e 2 A0 + /(0) +eM ](t'X (17) 

and 

V (t- T ) = t e~ r ( *-*' ) Au(t') dt, (18) 



where, again, it is understood that x an d v are also functions of the space variable x and 
are defined for all x G Q, we obtain 



u t + e 2 (j) t = J °° «i(r) v(t; r)dr, 



v t (t; r) + r v(t; r) = An, 



e 2 & = / °° a(r) X (t; r) rfr, 



(19) 



Xt(t; r) + r x(t; r) = e 2 A0 + /(0) + e u, 
for all t G [0,oo). Note that the second and fourth equations satisfy v(0, r) = x(0, r) = 
for all x G £1 and for all r G [0, 00). 

For an exponentially decreasing kernel @, 5(r) = <5(r — 7). For a kernel which is a linear 
combination of kernels of type (f|), a(£) = Z)fc=i a fe e~ 7fct , for positive constants and 7^ 
(A; = 1,2, .. . , n), we have dt{r) = J2k=i a k^{ T ~ Ik)- I n the latter case we can define for 
k — 1, 2, ... ,n 



X(*5 7fc) = / e 
Jo 

and substitute into ([I]) obtaining 



"7fe (*-f) [ J 



[e 2 A0 + /(0) + e/i ](*')< 



(20) 



e 2 0t = X(*5 7fc)> 



(21) 



+ /(0) +e/i, 



for = 1, 2, . . . , n. It is clear that fl2~T| ) is a particular (discrete) case of (IB) and the simplest 
generalization of ([$]). This discussion can be easily adapted to system (|5|). 



The approach discussed up to now in this section allows us to consider kernels which are not 



necessarily exponentially decreasing. Examples are, for positive a, (3 and 7: a(t) = l/(jt+a), 
for which a(r) = 1/^ e~~ l or, more generally a(t) = l/(7t + oc) n (with n a positive integer), 
for which a(r) = 1/(7" (n - 1)!) r"" 1 e - *'; a (t) = /3/((7 t + a) 2 + (3 2 ) and a(t) = ( 7 t + 
a )/{{lt + a ) 2 + /^ 2 ) f° r which a(r) = I/7 2 e~~*sin(/3 r) and a(r) = I/7 2 e~~*cos(/3 r) 
respectively. A class of kernels not included in (|l^) or ( |14|) is a(t) = a e~ 7 *cos(/3 t), or 
equivalently, a(i) = a e Z1 *+ e Zlt where Z\ — 7 — % (3. 

From a slightly different point of view, kernels of type a(t) = a e~ 7 * cos(/5 t) can be repre- 
sented as Fourier transforms of suitable functions. 



or 



00 

a(t) = I a(w) e~ iwt dw, (22) 



Oi(t) = I di(w) e~ iwt dw. (23) 

for i = l,2 where 



a(w) = — / a{t) e lwt dt, 



2ti 



CO 



and 



Oi(w) = — / <n{t) e lwt dt. 

-oo 



2tt 

for i = 1,2. In this case the functions a(t) and ctj(t), i — 1,2 are understood to be defined 
for all i and being equal to zero for t < 0. Substituting ( p2[) into ([I]), rearranging terms and 
defining 

X (t;w)= [ e 2 A + /(0) + e fc](f)df, (24) 

where it is understood that x is a function of the space variable 1 6 fi, we obtain 
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e 2 <k = I-oc a ( w ) X(t; w) dw, 



(25) 



_ Xt(t; w)+iw x(t] w) = e 2 A0 + f(<f>) + e h, 
for all w G (—00, 00). For the phase field equations with memory, substituting (|23|) into (|5|) 
rearranging terms and calling 



X (t; w) = / e"™^') [ e 2 A0 + /(0) + e u ](f) df 
Jo 



(26) 



and 



Jo 



(27) 



where, again, it is understood that \ and v are functions of the space variable xefl. System 
(|) becomes, for all w 



u t + e 2 (f) t = ai(w) v(t; w)dw, 



v t (t; w) + i w v(t; w) = Au, 



e 2 <Pt = I-oc a ( w ) X(t;w) dw, 



(28) 



Xt(t; w) + iwchi(t; w) = e 2 A0 + f(<p) + e u, 
for all w G (—00,0x3). As we set at the beginning of the introduction, the kernels a(t), are 
assumed to be piecewise continuous, differentiable at the origin, scalar-valued functions on 
(0, 00). Moreover they are assumed to be independent of e and such that 



/ a(t) dt < 00, 
Jo 



and 
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Jo 



■oo 



<y(s) ds < oo 



and / a(s) s ds 



< oo. 



Jo 



when we use the Laplace transform, and 



'OO 



a(s) ds < oo 



and / a(s) s ds 



< oo. 



— oo 



J —oo 



when we use the Fourier transform. In both cases this is equivalent to a(0) < oo and 
a'(0) < oo. 

3 The Allen-Cahn equation with memory - asymptotic 
analysis in Cartesian coordinates 

For points outside the interface and for e > we expand <fi as follows 



giving f'((f>°) (p 1 + h = 0. Thus, for points which have not yet been reached by the moving 
front, we have 



= 4>(x, t] e) = 4>°(x, t) + t t) + e 2 2 (x, t) + C(e 3 ), 



and substitute into (P obtaining the 0(1) and 0(e) respectively 




whose solution is f(<f>°) = (or 0° = ± 1), and 




(29) 



and 



12 



X \t;r) = 0, (30) 

for all r G [0, oo). 

For the asymptotic analysis using Cartesian coordinates the interface is represented by 

y = S(x, t, e) for e sufficiently small and assume St 7^ for all x G VL and all £ > 0. We define 

a new variable 

_ = |/ - 5(a;,t,e) 



which is assumed to be 0(1) as e — > near the interface. We call $ the asymptotic form of 
<p as e — ► with z fixed; i.e., 

(j> = $(z,x,t;e). (31) 

The field equations (|T6| ) in (z, x, £) coordinates, after differentiating the first equation with 
respect to t and calling \ = Y(x, z, £; r, e), become (see Appendix A) 

e 3 $ tt - 2 e 2 Si + e Sf - e 2 S tt = / °° a(r) [ e T* - S t T 2 ]dr, 

eT i -S t T 2 + f rT = f 3 $ n -2e 2 S, $ 2:E + e (1 + S 2 ) § zz - e 2 S xx $ z + e /($) + e fc, 

(32) 

for r G [0, oo). 

The asymptotic expansion of $ is assumed to have the form 

$~$ + e$ 1 + e 2 $ 2 + C(e 3 ), as e -> 0. (33) 

Substituting into (|32]) and equating coefficients of the corresponding powers of e, we obtain 
the following problems for 0(1), 0(e) and 0(e 2 ) respectively: 
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' J °° a(r) S t T° dr = 0, 

< 

k S t T° = 0, 
' S?*2. = Jo° «(r) [T?-5 t Ti dr], 

< 

_ T°-^Ti + rT° = (l + 5 a 2 ) + 

and 

' -2 5, $° t + 5? $L - Sit $ S = Jo 00 S(r) [ T, 1 - S t T\ ] dr, 

< 

T\ - S t Tl + r T 1 = -2 & $L + (1 + 52) *L " S m $° + f(<&°) $ 1 + h, 
for all r G [0, oo). From ( p9|) and fl30l) we have 

lim T°(x, z, t; r) = 0, 

z — »oo 

and 

lim T 1 (x, z, t: r) = 0. 

From fl54"|) we have T°(x, 2, i;r) = for all r G [0, oo). Integrating with respect to 
applying condition ( |3^ ) (assuming St ^ 0) yields 

T° = 0, 

for all r G [0, oo). Substituting fl39|) into fl55|) the second equation becomes 

-5 t Ti = (l + S2) + 
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We see that T\ does not depend on r. We multiply the second equation in ( p5| ) by a{j) and 
integrate with respect to r. We obtain 

POO 

^$L = [(l + ^)$L + /(«f°)] / a(r)dr 

Jo 

since <3> and S do not depend on r. We call 

a = g(r) d-r) 1 , (41) 

obtaining 

( 1 + Si - a S 2 t )$L + /($°) = 0. (42) 

From (f4"0|) we have that St T* = —a S| Integrating with respect to z and applying 
condition fl38T) we obtain 



T 1 = -a S t $°. (43) 
In order to solve ((421) we define a new variable 



£ ■= ~ r- (44) 

(1 + 5* -a 5?)* 



In terms of £, equation ((42]) reads 

d>° ? + /(d>°) = 0, (45) 

whose solution is $° = ), the unique solution of = 0, *(±oo) = ± , ^(0) = 0. 

Thus 



$o = $ o ? (46) 

^(l + S*-aSf)U 
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Multiplying the second equation in (^) by a(r), integrating with respect to r, replacing 
(^3|) into fl36|), and rearranging terms we obtain 



( 1 + Si - a 2 St W zz + /'($ u ) $ x = (S xx - a S U W, - 2 a S t ®» t + 

+ 2 S x & zx -h- a 2 S t (7°° a(r) r dr) (47) 

Setting 



7 = a 2 (J a(r) r rfrj . (48) 
In terms of £, x and t equation fl4"T|) reads (see Appendix A) 

$i + f'($0) $i = S ** ~ aS tt~ jSt $ o _ 2Sx(S x S xx -a S t S xt ) ^ 

^ V ' (l + S*-a 2 S?)l € (l + S*-a 2 St)l 1 K ^ 



+ — 77 1 * ; + *?) - (49) 

(1 + 52- «2 5t) 2 

It is straightforward to check that ^'(O satisfies the homogeneous equation 



+ /'($°) $ 1 = 0. (50) 
This implies that the operator A defined by 

A:=| + W0) (51) 
has a simple eigenvalue at the origin with as the corresponding eigenfunction. Now the 
solvability condition for equation fl4~9|) gives 



S xx - a 5tt - 7 5 4 /- 00 ATr/ 

(l + 52-Q! 2 ^ 2 



, /-oo /-oo 

r / (tf') 2 ^ - h / 'd£+ 

2 J —OO J — oo 
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2 S X (S X S xx — a St S x t) — 2 a St (S x S x t — a St Stt, 
(1 + SI - a 2 S]jl 

A simple calculation shows that 



/oo 
(f + *')*'^ = 0. (52) 
-OO 



/oo 1 roo 

&9»dt = -- (V) 2 d£. (53) 
-oo / J —oo 



Defining 



, ^(+oo) - *(-oo) /c 

* : ^ ft • < 54 > 

Substituting (|53|) and (|5^) into fl52|) , multiplying equation ( (52|) by (1 + S% — a 2 Sf)* and 
rearranging terms we obtain (|ll]). Note that for /(</>) = (<p — <f) 3 )/2 (Ginzburg-Landau 
theory), = tanh | and v := 3 h, whereas for /(</») = sin0, = 4 tan _1 e^ — 7r and 

1/ := f /i. 

4 The phase field equations with memory - asymptotic 
analysis in local coordinates 

4.1 Assumptions and definitions 

For the asymptotic analysis, using local coordinates, we set x = (xi,x 2 ) and call d(x) the 
distance from x to V; i.e., d(x) = dist(x,T). We next define a local orthogonal coordinate 
system (r, s) in a neighborhood of T in the following way 



r(x, t;e) — < 



d(x) if <j)(x) > 

(55) 

-d(x) if (j)(x) < 0, 
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and s(x,t;e), a smooth function of t, such that on T(t; e) it measures arclength from some 
point which moves normally to T as t varies. The assumed initial smoothness of F(t; e) 
implies that r is a smooth function, at least, in a sufficiently small neighborhood of T. 

The outer expansions of <fi and u are assumed to have the form 

= <j>(x, t; e) = (f)°(x, t) + e <f>\x, t) + e 2 2 (a;, t) + C(e 3 ) (56) 

and 

u = u(x, t; e) = u°(x, t) + e u x (x, t) + e 2 w 2 (a;, t) + 0(e 3 ). (57) 
In order to determine the inner expansions we first define the inner variable 

z(x,t;e) := T -^M ( 58 ) 
and then assume the inner expansions to be given by 

(f) = $(z, s, t; e) = $°(z, s,t)+e ^\z, s, t) + C(e 2 ) (59) 

and 

u = U(z, s, t; e) = U°(z, s,t) + e U\z, s, t) + C(e 2 ). (60) 

The very definition of V requires $(0, s,t; e) = 0. In what follows we will use the following 
notation to refer to any variable g evaluated by approaching T from either side (r > or 
r < 0): 

g | r ±= lim g(r, s,t;e), (61) 
g r | r ±= Mm g r (r,s,t;e). (62) 
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The following relations between the inner and outer variables obtained in [?] are assumed to 
hold as p — > ±00. 



G°(p,s,t)=g°(0 ± ,s,t), 



(63) 



G\p, s,t)=g\0 ± ,s,t)+pg° r (0 ± ,s,t). 



(64) 



4.2 Derivation of the equations of motion 

Let us look at the system fl5|) assuming that u = initially and along the boundary. 



4.2.1 Outer problems 



Substituting fl5"5|) and ( |5"TD into (|5|) and equating coefficients of the corresponding powers of 
e we obtain the 0(1) and 0(e) outer problems respectively for points where the interface 
has not yet arrived: 



and 



it? = ai * Au c 



a 2 * f((f>°) = 0, 



(65) 



u\ = a\ * Am 1 , 



(66) 



The solution of (p5|), given the assumed initially and boundary conditions for u, is m° = 0, 
6° = ±1. The solution of (1661) is w 1 = 0, ^ = 0. 
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4.2.2 Inner problems 

From the solution of the outer problems fl65|) and (|66| ) and for points where the moving front 
has not yet arrived, we have 

x °(t;w) = 0, and v°{t;w)=0, (67) 

and 

x \t]w)=0. and v 1 {t;w)=0, (68) 
for all x G Q and w G (— oo, oo). 

System (|28f) expressed in the (z, s,i) coordinates (see appendix B), after differentiating 
the first and third equations with respect to t and calling v = V(z,s,t;w,e) and x = 
T(z, s, t; w, e), becomes 

r 2 U zz + e[2r t U zt + r tt U z } + e 2 [U tt + 2 s t U st + s 2 U ss + s tt U s + 2r t s t U zs + r 2 <$> zz ] = 

= I-oo «iH [e 2 V t + er t V z + e 2 s t V s ) dw + 0(e 3 ), 
< er t V z + e 2 [V 4 + s t V s + i w V] = U zz + e Ar C/ 2 + e 2 [ Z7„|Vs| + ?7 s As ], 

er 2 $ zz + e 2 [ 2 r t $ 2t + r« $ z + 2 r t s t ] = Ho a 2 H[eT t + r t T, + e s t T s ] c^ + 0( e 3 ), 
/l T, + e [T ! + Si T s +iwT] = e[^ + /($) ] + e 2 [ Ar $ z + [/ ] + G(e 3 ), 

(69) 



for w E (—oo, oo). From (|67j ) and (pal) we have 



20 



lim T°(x, z, t; w) = 0, 



and 



lim V°(x, z, t: w) 



and 



lim T 1 (x, z, t: w) = 0, and lim V 1 (x, z, t: w) = 0. 

Z— >OC> 2— >00 

Substituting (|59|) and ( |60|) into (l69| ) and equating coefficients of the corresponding 
of e we obtain the 0(1), 0(e) and 0(e 2 ) problems respectively. 

0(1): 



r 2 r/° = 

't u zz u ' 



eft = o, 



n r° = o, 



0(6): 



rt K° = Ul + Ar I/J, 
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' r 2 Ul + 2 r t U] t + T tt U] + E/g + 2 s t V% + s 2 C/° + s tt L/J + 2r tS( l£ 

+rl = IZo «iH[ + r, V? + ^ V s ° } dw, 
* n V} + V? + s t V s ° + i wV° = Ul + Ar + [/°|Vs| + U° s As, (74) 

r t 2 $L + 2 r t + r« $° + 2 r t * <f° s = Jl^ a a («;)[ T, 1 + r t T 2 + s t ] dw, 
k TJ + r t X 2 + a* TJ + z w T 1 = $^ + f($°) "I" 1 + Ar $° + [7°, 

for w G (— oo, oo). We call 

a.i = {^J oa(w) dw^j , (75) 

for i — 1, 2. A bounded solution of the first two equations in (|72| ) satisfying the matching con- 
ditions (|63| ) is f/ 1 = 0. From the third and fourth equations in ( |72|) we have T°(x, 2, £; w) — 
for u> G (—00, 00). Integrating with respect to z (assuming that r t 7^ 0) and applying condi- 
tion ( [70|) yields 

T° = 0, (76) 

for u> G (—00,00). Substituting E/ 1 = in the first two equations in ([73|), multiplying the 
second and fourth equations in (|73|) by <ii and o7 2 respectively, and integrating with respect 
to w yields 
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(a, - r?) Ul = 0, 

(77) 

(1 - a 2 r?) + = 0, 

Assuming that r\ ^ a±, the bounded solution of the first equation in (|77|) satisfying the 
matching conditions (|63|) is f/ 1 = 0. To solve the second equation in (|77|), we assume that 
rf 7^ «2 and we define the new variable 



1 - a 2 n 2 )2 



(7f 



In terms of (£, s,t), equation the second equation in (|77|) reads 



+ /(*°) = 0, (79) 

whose solution is $° = #(£), the unique solution of = 0, *(±oo) = ±1, ^(0) = 0. 

Thus 



n ^ , (80) 

V(l-a 2 rf)2) 

which satisfies (0). From the fourth equation in (|73|) we have r t = a 2 r 2 <3>° 2 . Integrating 
with respect to z and applying condition (|7l| ) we get 



T 1 = -a a r 2 S° 



(81) 



Substituting fl30|) into the third and fourth equations in (|74D, multiplying the fourth equation 
by a 2 and integrating with respect to w yields 



;i - «2 r 2 ) $L + /' ($°) $ 1 



(«2 r« + 72 «2 r t - Ar) $° + +2 r t $ 



o 



52) 



23 



where 



72 = i a\ I J a2(w) w dwj . (83) 
Equation fl82|) expressed in the (£, s,i) coordinate system reads (see appendix B) 

+ /'($°) $ 1 = 



(1 - a, r?)» « " (1 - a 2 rf)i 
It is straightforward to check that ^'(0 satisfies the homogeneous equation 

+ f'($°) < & 1 = 0. That means that the operator A := + has a simple 

eigenvalue at the origin with ^' as the corresponding eigenfuction. The solvability condition 
for the equation (151) now gives 



2 rf r tt 



,1 - «2 r£ 
?"« + 72 r t - Ar 



/oo 



oo 

/\2 



(tf')X = 0. (85) 



oo 



(1 - a 2 r?)a 

A simple calculation shows that £ty'ty"d£ = —\ f^^') 2 ^. Hence multiplying equation 
(|35|) by (1 — «2 r|)t and rearranging terms one obtains 



— — 2 + 72 r t = «. (86) 



1 - «2 rj 

Taking into consideration that on the interface Ar = k, the curvature of the interface, and 
that r t = —v, its normal velocity [?], equation ( p6|) becomes (|T2|). 
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5 Conclusions 



In this paper we showed that to the leading order and for a large class of kernels a, a± and 
02 under suitable assumptions on them, the law governing the evolution of interfaces for the 
integro-differential equation ([!]) is the same as for the differential equation (^j. It is easy to 
see that 72 is given by 

_ a'(0) 

For equation (|5|) the result is similar with 7 and a changed by 72 and a 2 . Our derivation is 
only valid for advancing fronts; i.e., non-oscillating fronts in the sense that there exist points 
on the interface with vanishing velocity. This restriction prevent us from analyzing system 
where those oscillations may be relevant. This assumption is not necessary when the kernels 
are exponentially decreasing [?, ?]. 

To solve the 0(e) problem we have assumed that |r t | 7^ ^fcTy («i being similar to 0:2). On 
the other hand, assuming that 1 — «2 v 2 > initially and that there is not change of sign in 
this direction during the evolution, a fact which is true for the circular case, we can easily 
see from equation flE^ ) that \r t \ < Therefore, if a\a2 > 1 the former assumption does 
not add further restrictions on the interfacial motion. 

Acknowledgement: The authors are indebted to A. Novick-Cohen for the formulation of the 
problem and to Ch. Charach for valuable discussions. 

A Cartesian coordinates 

From Cartesian to (z,x,t) coordinates: 
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To go from Cartesian to (z, x, t) coordinates we transform derivatives as follows 



<k = &t - - S t $ x , 
e 



<f> tt =$ tt --2S ( $ zt + \ S] <5> zz - - S tt $z, 
e e 2 e 



1 1 1 

d = <f> 9 5 $ H S" 2 $ 5 $ 

"rxx ^ XX ^ U X ^ ZX ~ 9 u x zz ^XX ^ Zl 

e e 2 e 



^yy £ 2 < ^ >2 - 2 ' 



4>xxt — ^xxt S t &xxz 2 S xt & zx 2 S x & xzt H 9 2S X 5 t 
e e e e 2 



+— 2 S x S x t &zz H — k S x &zzt o S St &zzz S xxx $ 2 S xx &zt H — 9 S xx St 

e 2 e 2 e 6 e e e 2 



= _ $ £ $ 

9 ^ZZt q ^ ZZZ' 

e 2 e d 
From (z,x,t) to (£,x,t) coordinates: 

To go from (z, x, t) to (£, t) coordinates, derivatives are transformed as follows 

iz = (1 + S 2 X - a 2 S?)-3, 

^ = + 5 2 _ ^ ^-f _ ^ ^ 

& = -z(l + Sl-a 2 S^ 2 )" (S x S xt - aS t S tt ), 
izx = -(1 + S 2 X - a 2 St)~i (S X S XX - aS t S xt ), 
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U = -(1 + Si - a 2 S 2 t )~ {S x S xt - aS t S tt ). 



$o = = $o ( 87 ) 

2 5 (l + S*-a 2 S?)$ ? 1 ' 



*°* = + = - ,?f* aStS *s + (88) 



C = + *?U = - 7^7% + *?)■ (89) 

[i + b x - a 2 b t )2 



B Local coordinates 

B.l From Cartesian to (r,s,t) and (z,s,t) coordinates 

To go from Cartesian to (r, s, t) coordinates we transform derivatives as follows ( using the 
fact that |Vr| = 1 ) 

<j>t = $ t + $ r rt + ® s s t , 
(fa = $tt + 2$ rt r t + 2$ st s t + $ rr rf + $ ss s t 2 + $ r r tt + $ s s tt + 2$ rs r t s t , 

and 

A0 = $ rr + $ r Ar + $ ss |Vs| 2 + $ s As, 

where the operators V and A refers only to the spatial variable x. These derivatives expressed 
in terms of (z, s, t) are 
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<f) t = $ t + -®zTt + ®sS t , 

e 



(p tt = $ ti + 2-<S> zt r t + 2® st s t + -r$ zz r t 2 + $ ss s 2 + -$ 2 r tt + $ s s t t + 2-$ 2S rts t , 
e e z e e 

and 

11 -9 

A0 = + -$ 2 Ar + |Vsr + * 4 As, 

e 

B.2 From (z,s,t) to (£,s,t) coordinates 

To go from (z, s,t) to (£, s,t) coordinates, derivatives are transformed as follows 

& = (l-ar t 2 )-*, 

& = ^(1 - a r 2 )~5 artrt( 
£rt = (1 - a r 2 )"i ar t r tt . 



$° = = - 1 (90) 
(1-a rf)2 



* * = + = TT^V + (91) 

(1-a rf)2 
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